Adiabatic evolution under quantum control 
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One of the difficulties in adiabatic quantum computation is the limit on the computation time. 
Here we propose two schemes to speed-up the adiabatic evolution. To apply this controlled adia- 
batic evolution to adiabatic quantum computation, we design one of the schemes without any prior 
knowledge of the instantaneous eigenstates of the final Hamiltonian. Whereas in another scheme, 
the control is constructed with the instantaneous eigenstate that is the target state of the control. 
As an illustration, we study a two-level system driven by a time-dependent magnetic field under the 
control. The physics behind the control scheme is explained. 
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Adiabatic quantum computation (AQC) was first pro- 
posed in 2000 [lj as a means to solve NP-complcte prob- 
lems. It is polynomially equivalent to conventional (gate 
model) quantum computation Q and possesses some de- 
gree of fault tolerance. AQC can be described by the 
following Hamiltonian, 



H Q (t) = [l-X(t)]Hi + X(t)H f , 



(1) 



where the quantum system govern by Ho(t) evolves 
slowly with time t and remains in its ground state as 
X(t) changes monotonically from to 1 within a time 
T. The initial Hamiltonian Hi is assumed to have an 
easily accessible ground state into which the system is 
initialized, while the ground state of the final Hamilto- 
nian Hf encodes a problem's solution. In order to reach 
the final ground state with high fidelity, the adiabatic 
theorem requires T oc (AE min )- s , where AE 

min IS the 

minimum gap between the two lowest energy eigenstates 
of Hq{£). The power 5 can be 1, 2, or possibly some 
other number depending on the functional form of X(t) 
and the distribution of the higher energy levels [H, H, |j| • 
This leads to a limit on the computation time T, which 
holds true for adiabatic quantum computation without 
quantum control, thus it is natural to ask whether one 
can use quantum control to sjaeed up_ AQC? 
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about the application of classical and modern control 
strategy to quantum systems. It has generated increas- 
ing interest in the last few yea rs due to its potential En; 
plications in metrology 11 
and other technologies 



as well as 

its theoretical interest on its own. Several approaches 
to the control of a quantum system have been proposed 
in the past decade, which can be divided into coherent 
(unitary) and incoherent (non-unitary) control, accord- 
ing to how the controls enter the dynamics. In the co- 
herent control scheme, the controls enter the dynamics 
through the system Hamiltonian. It affects the time evo- 
lution of the system state, but not its spectrum, i.e., 
the eigenvalues of the target density matrix pf remain 



unchanged in the dynamics. In the incoherent control 
scheme, an auxiliary system, called probe, is introduced 
to manipulate the target system through their mutual 
interaction. This incoherent control scheme is of rel- 
evance whenever the system dynamics can not be di- 
rectly accessed and provides a non-unitary evolution for 
the quantum system. This breaks the limitation for the 
coherent control mentioned above. Among these quan- 
tum control strategies, Lyapunov control plays an im- 
portant role in quantum control theory. Lyapunov func- 
tions, originally used in control to analyze the stability 
of the control system, have formed the basis for new con- 
trol design. Several papers have be published recently to 
discuss the application of Lya pun ov control to quantum 
systems 0, [M, El E3, [H, S M EH. Although the ba- 
sic mathematical formulism for Lyapunov control is well 
established, many questions remain when one considers 
its applications in quantum information processing, for 
instance, whether one can use quantum Lyapunov con- 
trol to improve the adiabatic evolution, and consequently 
minimize the computation time? 

In this paper, we shall address this issue by using a two- 
level model with Lyapunov control. Since the two lowest 
levels are important for AQC, this model to some extent 
can good quantify the AQC under control. Two control 
schemes are proposed which correspond to two different 
choices of Lyapunov function. By numerical simulation, 
we find that these schemes work well. The paper is orga- 
nized as follows. In Sec. II, we present a general formalism 
for the Lyapunov control, two Lyapunov functions which 
will give two control schemes are constructed, then we use 
these schemes to manipulate a two-level system in Sec. III. 
Conclusion and discussions are presented in Sec. IV. 

General formalism. — Let us start by focusing on the 
Hamiltonian H$(t) in Eq.([T]). We denote the instanta- 
neous eigenstates and the corresponding eigenvalues of 
H (t) by \n(t)) and E n {t), respectively. AQC is designed 
to take advantage of the adiabatic theorem, it works by 
evolving a system from the accessible ground state of 
an initial Hamiltonian Hi to the ground state of a final 
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Hamiltonian Hf. For the system to remain in its ground 
state, the adiabatic theorem requires that the system 
Hamiltonian changes very slowly and there is an energy 
gap between the ground state and the others. The goal 
of this paper is to manipulate the system such that it 
remains in one of its eigenstates. To meet the require- 
ment of AQC, any information about the instantaneous 
eigenstates of Hf is forbidden to use in the control, but 
a prior knowledge regarding the system Hamiltonian is 
allowed. This motivates us to propose the first scheme 
in the following. Beside the interests in AQC, controlling 
a quantum system to a target state is interesting on its 
own. As the target state is known, the information about 
the target state can be used in the control design. This 
is the goal of the second scheme discussed in this paper. 

Scheme A. — We aim to design a control scheme that 
can manipulate a system to remain in its instantaneous 
ground state without any prior knowledge of its instan- 
taneous eigenstates of Hf. To this end, we introduce 
control operators H c j (j — 1,2,3,...) and require that 
[Hq, H C j] = for any j. The control operators H C j enter 
the system through control field f 3 {t). The total Hamil- 
tonian of the system is then written, 



H = Hn 



(2) 



where the control field fj(t) can be established by Lya- 
punov control theory. Define a function V(t) as 



v(t) = \m)\x\m)\ 2 



(3) 



with X being a time-independent Hermitian operator, we 
find V(t) > and 

v(t) = j{mm(t)Km\[H,x]m))- (4) 

Here and hereafter \ip(t)) denotes the state of the system 
that satisfy ih-^\^p(t)) — H\ip(t)). We now show how to 
establish the control field fj(t). Lyapunov control theory 
tells that for V(t) to be a Lyapunov function, V(t) has 
to satisfy, V > and V < 0. So, if we choose the control 
field fj(t) as 

= t(m\x\m)(m\[H CJ ,x]m)), for j+k, 

JjoW = - i ,i.\\\u FTTTTTu' for 3= Jo, 5 

(ip(t)\[H cjo ,X\\i){t)) 

then V < 0. Here jo was specified to satisfy 
{tp\[H C j , X)\tp) / 0. Clearly fj(t) is a time-dependent 
real number, thus the total Hamiltonian H is Hermitian. 
The key point of this control scheme is the choice of op- 
erator X, it dominates the success rate of the control and 
show the merit in this scheme. From the control design 
of the scheme, we find that the instantaneous eigenstates 
do not enter the control, this mean that by the present 



control strategy, we can manipulate the quantum sys- 
tem to evolve along one of the eigenstates without any 
knowledge of its instantaneous eigenstates. This is ex- 
actly what we want in AQC. To see that this control 
scheme indeed works, we note that [Hq,H c j] — 0, indi- 
cating the control H C j itself can not induce population 
transfer among the instantaneous eigenstates. Suppose 
that the system is initially prepared in its ground state 
|0(t = 0)) and H cj = H a , Eq.© yields f jo (t = 0) = -1 
and fj(t = 0) = 0, (j = 1,2,3, j f j ), regardless of 
what X takes. When \ip(t)) deviates from \0(t)), fj(t) 
can be very large depending on operator X. The Lya- 
punov control will then render the system nonlinear, and 
this nonlinear effect would bring the state \^(t)) back 
to |0(i)). We note that -gf- (as well as can derive 

the system from one instantaneous eigenstate to the oth- 
ers, this together with the control keep the system in the 
instantaneous eigenstate into which the system was ini- 
tially prepared. We will demonstrate this point through 
an example in detail later. 

Scheme B. — For general control problem, the target 
state is known, we then can use the target state to design 
a Lyapunov function and establish the control field fj (t) . 
Suppose that the target state is the zeroth instantaneous 
eigenstate \0(t)) of Ho(t), define 



v(\o(t)),m))) = i-\(o(t)m))f 



(6) 



Clearly V > with equality only if |0(i)> = IV'C*))* To see 
Eq.® indeed defines a Lyapunov function, we calculate 
the time derivative of V a,s(H = Ho 
scheme) , 



^ (i)^. in this 



v = -^E^W(^( < )l^|oW)<o(t)|^(t))-c.c.) 

3 

- 2Re«6(t)|^)>(V(*)|0(*)». (7) 
Obviously, 

fj(t) = -2Im(m)\H' cj \0(t)){0(t)\i,(t)}), for j ? j , 

f M _ hM(b(t)W))m)W))) _ . m 

/ioW Im((0(t)|^, o |^(t))(^|0(t)))' 3 - J0 ^ } 

guarantee that V < 0. Again, j n was selected to satisfy 
Im((0(t)|F^- o |V , (*))(V'|0(t))) ^ 0. Hence, the evolution of 
the system with Lyapunov control can be described by 
the following nonlinear autonomous equations, 

in^(t)) = H\m), 

H = H +J2f j (t)H' cj , 

3 

H o \0(t))=Eo\0(t)), (9) 

where f (t) (j = 1,2,3,...) are given by Eq.@. Note 
that [H' C j,Ho] = is not required in this scheme. The 
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difference between the present scheme and those in the 
literature [3(| is the target state. In earlier studies, the 
target state \<pj) is either time-independent or time- 
dependent via iU-§f\4>d) — Hq\4>^)^ whereas in our scheme, 
the target state is one of the instantaneous eigenstate 
of Ho. The time derivative of the instantaneous eigen- 
state plays an important role in the control, leading to 
totally different control fields fj(t) in Eq.®. Note that 
the choice of fj(t) in Eqs.JU) and © are not unique. 
In fact, when there is only one control operator H c , the 
control field f(t) in Eq.© can be chosen as, 

/(*) = i(m\x\m)(m\[H c ,x]m)) 
_ m)\[H ,x]m)) nm 

provided f(t) is a finite number. This is exactly the case 
in the example that we will illustrate below. 

Example. — As an illustration of the Lyapunov control 
scheme, we discuss below a two-level system driven by 
a time-dependent magnetic field. The Hamiltonian that 
describes such a system can be written as 

H (t) = fiB-a (11) 
= fiB (sin 9 cos <f>a x + sin 9 sin (f>a y + cos 9a z ) , 

where a x ,a y and a z are the Pauli matrices, Bq is the 
amplitude of the classical field. 9 is specified to be 
a constant here, while (f> depends on time t through 
<j> = ut (u>, constant). In comparison with Eq.|T]), Hi 
takes fiBo(sm9a x + cos 9a z ), while Hf = /j,Bo(sm9a y + 
cos 9a z ). Although in this model we can not found an 
analytical X(t) by which we write Eq. ljlip in the form of 
Eq.([T]), the model Eq. (fTTj) can be obviously mapped into 
Eq. ([!} . It is well know that the instantaneous eigenstates 
and the corresponding eigenvalues of Ho are \E+(t)) = 
cosfe^l T)+sinf| |), |£ (t)) =-sin|e^| T)+cos|| |) 
and E± — ±/ii?o, respectively. In the absence of Lya- 
punov control, it is required that lu -C [iB® for the system 
to evolve adiabatically. We now show that the system 
can evolve along one of its instantaneous eigenstates, e.g. 
|i?_(t)), under the Lyapunov control even if u) > fiBo. In 
the following, a fidelity defined by 

F(t) = \{E.{t)\m)\ 2 , (12) 

will be used to measure the effectiveness of the control. 
The dependence of F(t) on time t and the control oper- 
ator H' c as well as X are shown. The results show that 
these schemes work good with properly chosen X and 

K 

We first consider the scheme A, where no informa- 
tion about instantaneous eigenstates enter the control. 
Choose X = a x + Ra z with a rate R, the dynamics of 
the two-level system under control is governed by 

ih^-Mt)) = mm), 
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FIG. 1: (Color online) Fidelity defined by F(t) = 
\(E-(t)\ip(t)}\ 2 as a function of time. The system is initially 
prepared in \E-(t = 0)}. (a)-(e) correspond to different ratio 
R. From (a) to (e), R = 12,9,6,3,0, respectively. By con- 
trast, the result without control is shown by the dashed line. 
We have set /j,Bo — 1 for this plot. The time is shown in units 
of l/(/ii?o), the other parameters chosen are uo — 4 (in units 
of fiBo) and 9 = ?. Note ui = 4/i£?o, the adiabaticity breaks 
for the system without control. 

H = vB-a + f(t)H Cl 

fit) = i(m\xm))(m\[H c ,x]\m) 

- vBo. (13) 

Here H c was specified to be H c = (sin 6 cos <pa x + 
sin 6 sin (f>a y + cos 9a z ) in order to satisfy [H C ,H ] = 0. 
We have perform extensive numerical simulations for the 
dynamics Eq. (fl"3)) . selected results are plotted in Fig. Q] 
Two observations can be found from FigQ] (1) The 
scheme works only for special X. For some choices of 
X, the control favors the adiabatic evolution, but for the 
other X, the control spoils the adiabaticity of the evolu- 
tion. (2) For the specified X, X = a x + Ra z , the larger 
the ratio R is, the better the fidelity for the system in 
|£L(*)>. 

In the scheme B, a prior knowledge of the instanta- 
neous eigenstates is known and allowed to use to design 
the control field f(t). Choose \E-(t)) as the target state, 
the dynamics of the system is governed by Eq.© with 
replacing |0(t)) by \E-(t)). We show in Fig[^]the fidelity 
as a function of time with H' c = a z + ra x . Similar to the 
results given in the scheme A, the choice of the control 
operator H' c dominates the effectiveness of the control. 
Not all control operator H' c can help the adiabatic evo- 
lution. The difference between these two scheme is that 
the scheme B can reach fidelity 1 at the final stage, this 
may depend on the example demonstrated. 

For linear quantum system, the population transfer 
from one instantaneous eigenstate to the others depends 
on the energy gap and the operator that induces the 
transfer. The energy gap separates the instantaneous 
eigenstates and ensures the adiabaticity of the evolution. 
The larger the energy gap is, the higher the probabil- 
ity with which the system remains in the instantaneous 
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Time 

FIG. 2: (Color online) Fidelity versus time t given by the 
scheme B. The fidelity was defined as the same as in Fig[T] 
From (a) to (e), the ratio r takes 0,2,4,6,8, respectively. 
Dashed line shows the result without control. All quantities 
are chosen and plotted in the same units as in Fig[T] 
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FIG. 3: (Color online) Eigenvalues of Hamiltonian H as a 
function of time. The parameters chosen are the same as 
in Fig[Tj.(a), indicating that the system remains well in one of 
the instantaneous eigenstates under control. The inset shows 
an enlarged level anti-crossing point. 



eigenstate. In the scheme A, the operator that induces 
the population transfer is SH = + /(i)^), small 
8H would benefits the adiabaticity of the system. For 
nonlinear system, however, the nonlinearity plays an im- 
portant role in addition to the factors for linear system 
mentioned above. All these factors together lead to the 
observed features and can be used to understand the 
physics behind the features. 

To make the explanation clear, we plot the eigenval- 
ues for the total Hamiltonian in scheme A as a function 
of time in Fig|3l The same parameters as in Fig[T]-(a) 
are chosen, which ensure that the state of the system 
remains in the instantaneous eigenstate |£L(t)). Recall 
that the energy gap between two instantaneous eigen- 
states is 2 (in units of fiBo) for H , we find from Fig[3] 
that the energy gap is intensively enlarged at most times 
by the control. Nevertheless, the enlarged energy gap 
itself can not explain why the system remains well in 
the instantaneous eigenstate, because the energy gap is 



very small at the level anti-crossing points, this indi- 
cates that the nonlinearity must play an important role 
in the dynamics. Indeed our numerical calculation shows 
that at these point, the nonlinear coefficient defined by 
\(E-(i)\f(t)H c \E-(t)) \ is much larger than the tunneling 
coefficient \(E-(t)\6H\E + (t))\, leading to the observed 
feature reminiscent of self-trapping in nonlinear system. 

To sum up, we have proposed two schemes to speed- 
up the adiabatic evolution. In the scheme A, the control 
has been designed without any information about the in- 
stantaneous eigenstates of the final Hamiltonian, hence 
this scheme can be used in the adiabatic quantum com- 
putation. The scheme B has been proposed using the in- 
stantaneous eigenstate of the system Hamiltonian. This 
scheme is applicable to control a quantum system when 
the target state is known. To show how the schemes 
work, we have presented an example consisting of a two- 
level system in a rotating magnetic field. The fidelity to 
quantify the effectiveness of the scheme was calculated 
and discussed. The results show that the fidelity depends 
sharply on the choice of X and the control operator H' c . 
The physics behind the schemes is revealed, which can be 
understood as an effect combining nonlinear effects and 
the broadened energy gap. 
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